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Objectives: This study aimed to define the molecular basis of co-trimoxazole resistance in Malawian pneumococci 
under the dual selective pressure of widespread co-trimoxazole and sulfadoxine/pyrimethamine use. 

Methods: We measured the trimethoprim and sulfamethoxazole MICs and analysed folA and folP nucleotide and 
translated amino acid sequences for 143 pneumococci isolated from carriage and invasive disease in Malawi 
(2002-08). 

Results: Pneumococci were highly resistant to both trimethoprim and sulfamethoxazole (96%, 137/143). Sulfa- 
methoxazole-resistant isolates showed a 3 or 6 bp insertion in the sulphonamide-binding site of folP. The trimetho- 
prim-resistant isolates fell into three genotypic groups based on dihydrofolate reductase (encoded by folA) 
mutations: Ile-100-Leu (10%), the Ile-100-Leu substitution together with a residue 92 substitution (56%) and 
those with a novel uncharacterized resistance genotype (34%). The nucleotide sequence divergence and dN/dS 
of folA and folP remained stable from 2004 onwards. 

Conclusions: S. pneumoniae exhibit almost universal co-trimoxazole resistance in vitro and in silico that we believe 
is driven by extensive co-trimoxazole and sulfadoxine/pyrimethamine use. More than one-third of pneumococci 
employ a novel mechanism of co-trimoxazole resistance. Resistance has now reached a point of stabilizing evolu- 
tion. The use of co-trimoxazole to prevent pneumococcal infection in HIV/AIDS patients in sub-Saharan Africa 
should be re-evaluated. 
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Introduction 

In sub-Saharan Africa, Streptococcus pneumoniae is one of the most 
common causes of meningitis and pneumonia. The combination of 
trimethoprim and sulfamethoxazole (co-trimoxazole) is recom- 
mended by WHO for prophylaxis in HIV/AIDS patients to prevent 
opportunistic bacterial infections and Pneumocystis jirovecii. 1 
Co-trimoxazole preventative therapy (CPT) became national policy 
in Malawi in 2005 and is now administered across all government- 
funded hospitals. In 2006, 34942 Malawians were registered to 
receive CPT; by the end of 2010 this figure had increased to 
189 520. 2 Between 1993 and 2007, sulfadoxine/pyrimethamine 
was the recommended first-line treatment for uncomplicated 



malaria in Malawi and was used to prevent malaria in pregnancy. 3 
Prior to the introduction of CPT, co-trimoxazole resistance in 
S. pneumoniae in Malawi was already high at 74% and since 2005 
resistance has been consistently above 90%. 4 

Co-trimoxazole and sulfadoxine/pyrimethamine both target 
dihydrofolate reductase (DHFR) and dihydropteroate synthase 
(DHPS), allowing cross-resistance mechanisms to these two 
antimicrobials. 5 DHFR and DHPS form part of the folic acid biosyn- 
thetic pathway. Co-trimoxazole and sulfadoxine/pyrimethamine 
act as false substrate inhibitors, preventing folic acid biosynthesis 
and bacterial cell growth. Resistance to co-trimoxazole and sulfa- 
doxine/pyrimethamine is conferred by acquisition of mutations 
in folA and folP, the genes encoding DHFR and DHPS, respectively. 
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Although multiple mutations have been reported in trimethoprim- 
resistant pneumococci, the single substitution Ile-100-Leu in DHFR 
is sufficient to confer resistance. 6,7 In sulfamethoxazole- 
resistant isolates, folP is characterized by a 3 or 6 bp insertion, 
resulting in the insertion of one or two amino acids in the DHPS 
sulphonamide-binding site. 8 " 10 

We analysed the folA and folP sequences of 143 pneumococci 
isolated from carriage and invasive disease following the introduc- 
tion of CPTin Malawi, to define the molecular basis of co-trimoxazole 
resistance under the dual selective pressure of widespread co- 
trimoxazole and sulfadoxine/pyrimethamine use. 

Methods 

Ethics 

We performed a detailed characterization of bacterial isolates from clinical 
specimens taken from patients for clinical reasons. The isolates used in the 
study were anonymized. These data are published with the approval of the 
University of Malawi College of Medicine Research & Ethics Committee and 
conform to institutional guidelines. 



Study isolates 

The Malawi-Liverpool-Wellcome Trust Clinical Research Programme (MLW), 
based at Queen Elizabeth Central Hospital (QECH) in Blantyre, Malawi, has 
archived over 5000 pneumococcal isolates since 1996. A convenience 
sample of 143 pneumococci, collected between 2002 and 2008, was 
selected from the archive as part of a study into the genetic diversity of 
S. pneumoniae prior to the roll-out of 13-valent pneumococcal conjugate 
vaccine in Malawi. 11 The selected isolates encompassed 35 serotypes 
and 68 multilocus sequence typing (MLST) sequence types (STs); 12 65 iso- 
lates belonged to one of 48 MLST STs that had not previously been reported 
to the MLST database and are here referred to as STunknown, followed by a 
number. One hundred and thirty-one pneumococci were from invasive 
disease (92% of the total) and 1 2 were from nasopharyngeal carriage (8%). 



Culture and extraction methodologies 

Co-trimoxazole susceptibility was assessed using the BSAC disc diffusion 
method. 13 Sulfamethoxazole and trimethoprim MICs were determined by 
the Etest (AB Biodisk, Solna, Sweden) according to the manufacturer's 
instructions. S. pneumoniae ATCC 49619 was used as a quality control 
strain and gave values within an acceptable reported range. 

S. pneumoniae isolates that had been stored in Microbank bacterial pre- 
servers (Prolab Diagnostics, Ontario, Canada) were streaked onto blood 
agar plates and incubated at 37°C for 18 h. Single isolated colonies were 
suspended in Todd -Hewitt Broth (Oxoid, Basingstoke, UK) and incubated 
at 37°C for 18 h. The cells were sedimented by centrifugation and resus- 
pended in 480 |j± of 50 |j,M EDTA and 120 jjlL of lysozyme and incubated 
at 37 C for 1 h. Genomic DNA was prepared from the lysate using a 
Promega Wizard Genomic DNA purification kit (Promega, Madison, USA). 
Multiplex DNA sequencing was performed with an Illumina Genome Ana- 
lyzer GAII (Illumina, CA, USA), as described elsewhere. 14 All of the sequence 
reads generated are deposited in the short read archive (National Centre 
for Biotechnology Information) under the accession numbers ERP000185 
and ERP000152. Reads were assembled using Velvet vl.0.03, 15 and con- 
tiguated against a complete reference sequence (accession number 
FM211187) using ABACUS. 16 Serotype and MLST STs were determined as 
previously described. 17 The folA and folP genes were identified by BLAST 
searching fo!A_R6 (Gene ID 4442919) and fo!P_R6 (Gene ID 4443057) 
from the fully susceptible laboratory reference strain S. pneumoniae R6 (tri- 
methoprim MIC 2 mg/L; sulfamethoxazole MIC 16 mg/L). Significant hits 



Table 1. Amino acid variation of DHPS from sulfamethoxazole-resistant 
pneumococci in relation to the co-trimoxazole-susceptible laboratory 
reference S. pneumoniae R6 







Number of 


Previously reported 


folP insertion 


DHPS mutation 0 


isolates (%) 


in S. pneumoniae 


Wild-type 




5(3) 


30 


6 bp 


STRPRPGSSYVEIE 


88 (65) 


8,9,30 


3 bp 


STRPGR_SSYVEIE 


21 (15) 


8 


6 bp 


STRPSSSSSYVEIE 


15 (10) 


30 


6 bp 


STRPGSVSSYVEIE 


2(1) 


10 


6 bp 


STRPGTGSSYVEIE 


6(4) 




6 bp 


STRPGSSYVEVEIE 


2(1) 




6 bp 


STRPGSSYVEIEIE 


4(3) 


8,10,31 



Amino acid insertions are in bold italics. 

Represents amino acid residues 56 -67 in the amino acid sequence of DHPS 
from S. pneumoniae R6. 
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Figure 1. Diversity of foiP from Malawian pneumococci, 2003-08. 
(a) Annual sequence divergence of folP from pneumococci in relation to 
the co-trimoxazole-susceptible laboratory reference strain S. pneumoniae 
R6. (b) Graph showing the annual dN/dS of folP from pneumococci in 
relation to the co-trimoxazole-susceptible laboratory reference strain 
S. pneumoniae R6. 



369 



Cornicket oi 



— 1 100 



ft 




106_STu802 
107 ST802 
103 STu35 
104IST802 
134 ST802 
105 ST802 
132 STu33 
74_5Tu48 • 
143_STu47 

82 ST63 

79 ST63 
77IST63 
59 STU32 

80 ST63 

83 ST2678 
81_ST63 
78 ST63 
127 STu7 
96 5Tu2 
147_STu2 
88_STu3 
43 ST2987 
53 STu5 
67_ST989 
71_STul4 
156 STul6 
65_STul7 ■ 
70_ST989 

69 ST989 
160 ST989 
2_ST2927 
159 ST1047 
126 ST700 . 
94 5Tu22 * 
149_STul3 
116 STu9 
31_STul3 
20 STul3 
33_STul3 
32 STul3 
57 ST102 
92 ST705 
91_ST705 
STu56 
„j ST5902 
90 ST705 
R6" 

75 STU48 • 
51_STul7 ■ 



62_ST347 A 
170_STull 
123 STu20 
55 5Tul8 
44ISTu45 
58 STu24 
76 STU48 • 
54ISTu29 
119 STu44 

121 STU42 

122 STul9 
102_STu54 
117_STu55 

I STu8 
21_ST217 

25 STU23 
4JT217 

26 ST217 
131 ST303 
15JT217 
142_ST217 

157 ST217 

158 ST217 
19_ST217 
20_ST217 
28 ST217 
12 ST217 
13_ST217 

9 ST217 

10 ST217 
165_ST217 
8_ST217 

5 ST217 

6 ST217 
155_ST217 

17 ST217 

27 ST217 
14 ST217 
163_ST217 
16 ST217 
7JT217 

II ST217 
22_ST217 
3 ST217 
157 ST303 

18 ST217 
63_ST347 A 
169 STu31 

4100 4 5_STu*6 
' 97_STu57 
38_ST289 
41 ST289 

39 ST289 
163 ST289 
98_5Tu21 

40 ST289 
37 ST289 
35_ST289 
36_ST289 
24 ST289 
66ISTu6 
113 ST172 
112_ST172 
95 STu22 it 
141 ST2650 
146_ST2650 
115 STU26 
114_ST361 

135 ST361 
48_STu41 
47 ST2285 
49 STu34 
105 STu43 
42_5tu52 
109 ST2987 
168 ST1871 

136 STuSl 
140_ST2213 
93_ST2213 
133 STU12 
148 STul2 
29_5Tul2 
166 STU12 
150_STul2 
161 STul2 
52 ST2902 
85_STul 
144 STulO 
87 STul 
86_STul 
84 STu39 
100 STU40 



'_STu27 



were viewed, edited and annotated in Artemis V11.22. 18 Nucleotide 
sequences of folA and fo/P were aligned separately in Seaview V4.1 19 
using Muscle V3.7. 20 The nucleotide sequences were translated into 
amino acid sequences and a second alignment was performed. 

The ratio of non-synonymous to synonymous single-nucleotide poly- 
morphisms (dN/dS SNPs) in folA and fo/P for each isolate was computed 
using the Nei-Gojobori method. 21 

Phytogeny construction 

Phylogenetic trees were constructed using RAxML v.7.0.4 22 using the nu- 
cleotide alignments. A generalized time-reversible model with gamma cor- 
rection for among-site variation was used with 10 iterations using different 
starting trees. To assess support for relationships in the tree, 100 random 
bootstrap replicates were performed. folA_R6 or fo/P_R6 was included as 
an outgroup. 

Results 

In vitro resistance 

Co-trimoxazole resistance amongst the study isolates was 96% 
(137/143). The MICs of sulfamethoxazole and trimethoprim for all 
resistant isolates were >1024 and >32 mg/L, respectively. This is 
very high-level resistance to both antimicrobials based on BSAC defi- 
nitions. 13 Six isolates were susceptible to co-trimoxazole. Five co- 
trimoxazole-susceptible isolates (53, 88, 96, 127 and 147) were sus- 
ceptible to both trimethoprim and sulfamethoxazole (MICs 12 and 
96 mg/L, respectively). The remaining co-trimoxazole-susceptible 
isolate, isolate 137, was susceptible to trimethoprim but resistant 
to sulfamethoxazole (MICs 6 and >1024 mg/L, respectively). 

Molecular markers of sulfamethoxazole resistance 

All sulfamethoxazole-resistant isolates displayed an insertion of 3 
or 6 bp in fo/P, resulting in the insertion of one or two amino acids 
between Arg-58 and Ile-66 of DHPS, an area hypothesized to 
form the sulphonamide-binding site. 7 The five sulphonamide- 
susceptible isolates did not display amino acid insertions in this 
area (Table 1). No additional mutations were detected that were 
unique to the sulphonamide-resistant isolates. 

The mean nucleotide sequence divergence of fo/P, in relation to 
fo/P_R6, was relatively stable over the isolate collection period, 
from 5.66% in 2004 to 5.72% in 2008 (P=0.6) (Figure la). The 
mean dN/dS ratio of fo/P also remained stable over the same 
period, from 0.10 in 2004 to 0.12 in 2008 (P=0.09) (Figure lb). 

A maximum likelihood tree constructed using the nucleotide se- 
quence of fo/P from the 143 study isolates is presented in Figure 2. 
With the exception of isolate 127, all of the sulfamethoxazole- 
susceptible isolates formed a single clade separate from the 
trimethoprim-susceptible folA _RS. All of the serotype 1 isolates, re- 
gardless of their ST, shared an identical fo/Psequence and clustered 
in a distinct clade. Of the 20 STs represented more than once in the 
study population, 16 shared identical fo/P sequences within all 



Figure 2. Maximum likelihood phylogenetic tree based on the fo/P SNPs 
from 143 Malawian pneumococci, annotated with serotype and ST. STs 
represented more than twice in the study population that did not cluster 
on the same clade are highlighted by a symbol: STunknown22, star; 
ST347, triangle; Stu48, square; STunknownl7, circle. Bootstrap values are 
shown on each branch. The scale bar represents the number of SNPs. 
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isolates belonging to that ST; 16 separate clades were formed in 
the phylogeny, each consisting of isolates with only a single ST. 
The remaining four STs did not share identical fo/P sequences, 
and the isolates belonging to each ST sat on a distinct lineage 
within the phylogeny. Three serotype 7F, STunknown48 isolates 
were represented in the study collection (isolates 74, 75 and 76). 
Isolates 75 and 76 shared identical fo/Psequences and grouped to- 
gether in the same clade; however, isolate 74 sat in a distinct clade 
with ST16, ST63, ST6872 and STunknown32. The fo/P sequence of 
isolate 74 was 6.8% divergent from that of isolates 75 and 76. 
Isolate 74 had an fo/P sequence identical to that of the four 
other STs belonging to the same clade. 

Two of three serotype 12B STunknownl7 isolates, 50 and 51, 
grouped in the same clade; however, a third, isolate 65, sat in a 
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Figure 3. Diversity of fo/A from Malawian pneumococci 2003-08. (a) 
Annual sequence divergence of folA from Malawian pneumococci in 
relation to the co-trimoxazole-susceptible laboratory reference strain S. 
pneumoniae R6. (b) Graph showing the annual dN/dS of folA from 
Malawian pneumococci in relation to the co-trimoxazole-susceptible 
laboratory reference strain S. pneumoniae R6. 



clade with serotype 12B ST989. Isolates 51 and 50 had an identical 
fo/P sequence, but the sequence diverged from that of the fo/P of 
isolate 65 by 7.2%. The two 18B, STunknown22 isolates (isolates 
94 and 95) represented in the study population did not sit together 
on the tree. There was 5.6% sequence divergence between the fo/P 
from the two isolates. The two serotype 19F, ST347 isolates also sat 
apart in distinct clades (7.3% sequence divergence). The nucleo- 
tide sequence of fo/P from isolate 62 was unique within the study 
population, whilst isolate 63 had an identical fo/P sequence to 
isolate 12B, STunknownl7, which sat on the same clade as 
isolate 63. 



Molecular markers of trimethoprim resistance 

Nucleotide sequences of fo/A from all isolates were determined 
and compared with fo/A_R6. There was an annual increase in se- 
quence divergence from 2.82% in 2004 to 4.75% in 2008; 
however, this increase was not significant (P= 0.125) (Figure 3a). 
The dN/dS ratio of folA could not be calculated for 12 isolates as 
these isolates only displayed non-synonymous SNPs within the 
gene. The dN/dS ratio was calculated for the remaining 131 iso- 
lates. There was an annual decrease in the dN/dS ratio of fo/A, 
from 0.3 in 2004 to 0.2 in 2008; this decrease was not significant 
(P=0.72) (Figure 3b). 

Amino acid sequences of DHFR were predicted and analysed. 
Amino acid substitutions at two residues were dominant within 



Table 2. Amino acid variation at residues 92 and 100 of DHFR 
from trimethoprim-resistant pneumococci in relation to the 
co-trimoxazole-susceptible laboratory reference S. pneumoniae R6 



DHFR mutation 


Number of isolates (%) 


Asp-92-Arg, Ile-100-Leu 


54(41) 


Asp-92-Arg 


41 (30) 


Asp-92-Thr, Ile-100-Leu 


19 (14) 


Asp-92-Gly, Ile-100-Leu 


10 (8) 


Asp-92-Val, Ile-100-Leu 


8(6) 
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Figure 4. Annual distribution of three resistance genotypes amongst 137 
trimethoprim-resistant Malawian pneumococci. 
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the study population, I100L and substitutions at residue 92 
(Asp-92 to Ala, Gly, Val or Thr) (Table 2). The trimethoprim- 
resistant isolates could be assigned to one of three genotypic 
groups: those containing Ile-100-Leu (10%); those containing 
Ile-100-Leu in combination with a substitution at residue 92 
(56%); and those containing a substitution at residue 92 with 
the wild-type amino acid at residue position 100 (34%). In this 
third group, the amino acid substitution at residue 92 was 
always Asp-92-Ala. 

Other amino acid substitutions that were unique to the tri- 
methoprim binding site of the resistant isolates were Pro-70-Leu, 
Ala-78-Thr, Glu-94-Asp and Leu-135-Phe. Six trimethoprim- 
resistant isolates contained Asp-92-Ala, with no other amino 
acid substitutions in the trimethoprim-binding site. 

The distribution of three genotypic groups amongst the 
trimethoprim-resistant pneumococci on an annual basis is 
shown in Figure 4. A trend analysis comparing all of the study 
years showed a significant increase in the number of pneumococci 
possessing the dual genotype overthe study period (P= 0.045). The 
change in the number of pneumococci possessing the Ile-100-Leu 
phenotype was not significant (P=0.67). There was a decrease in 
the number of pneumococci with the Asp-92 mutations; 
however, this was not significant (P=0.06). 

A maximum likelihood tree constructed using the 138 variable 
sites in folA from the 143 study isolates is presented in Figure 5. 
The trimethoprim-susceptible isolates did not cluster together on 
the tree, sitting on four distinct clades. The ST was a poor indicator 
of how isolates would cluster on the tree. Only three STs (STun- 
known48, ST705 and STunknownl3) shared identical folA 
sequences and clustered together on the tree. The remaining 17 
STs did not cluster together and were dispersed across multiple 
fo/P lineages within the phylogeny. Isolates belonging to ST289 
were represented in four distinct folP lineages. 
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Discussion 

We propose that extensive use of both co-trimoxazole and sulfa- 
doxine/pyrimethamine in Malawi has been a major driving force 
behind the evolution and maintenance of co-trimoxazole resist- 
ance in Malawian pneumococci. Resistance amongst pneumo- 
cocci to co-trimoxazole was already high (74%) before the 
introduction of CPT in 2002. A previous study of Malawian children 
aged <5 years found that sulfadoxine/pyrimethamine treatment 
led to a significant increase in the proportion of children colonized 
with co-trimoxazole-non-susceptible pneumococci. 23 After the 
introduction of CPT, resistance rose to 92%. It has previously 
been reported that there is a significant correlation between re- 
gional consumption of co-trimoxazole and resistance in pneumo- 
cocci in the following year 24 and we believe the increased 
resistance is a result of extensive use of co-trimoxazole locally. 

All of the sulfamethoxazole-resistant isolates possessed a 3 or 
6 bp insertion in the nucleotide sequence encoding amino acids 



Figure 5. Maximum likelihood phylogenetic tree based on the folA SNPs 
from 143 Malawian pneumococci, annotated with serotype and ST. 
Trimethoprim-susceptible isolates are highlighted with a star. Bootstrap 
values are shown on each branch. The scale bar represents the number of 
SNPs. 
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58-67 of DHPS. This has previously been described as a hot spot 
for mutations in sulfamethoxazoie-resistant pneumococci. 3,9,10 
The most common insertions led to the duplication of Arg-58 
and Pro-59 and the insertion of an arginine residue between 
Gly-60 and Ser-61. These mutations have been previously shown 
to confirm resistance. 8 Two mutations, the Val 6 4Glu 65 and 
Gly 60 Thr 6 i repetitions described here, are novel. The role of these 
mutations in sulfamethoxazole resistance has yet to be estab- 
lished; however, their location in the sulphonamide-binding site 
amongst resistant isolates strongly suggests that these mutations 
reduce the affinity of sulfamethoxazole for DHPS. 

The phenotypic and genotypic characterization presented here 
shows that Malawian pneumococci are almost universally resistant 
to sulfamethoxazole. Previous studies have shown that insertions in 
the sulphonamide-binding site leading to a drastic reduction in the 
binding affinity to sulfamethoxazole have minimal impact on the 
binding affinity of DHPS to its natural substrate. 10 This suggests 
that, despite removal of sulfadoxine/pyrimethamine selective pres- 
sure on folP, sulfamethoxazoie-resistant pneumococci will continue 
to dominate within the population as the resistance imposes 
minimal fitness cost. 

Clonal dissemination of specific STs appears to be responsible 
for the dominance of sulphonamide resistance in Malawi. Phylo- 
genetic analysis demonstrated that ST was a strong indicator of 
isolate grouping. The 6.8% sequence divergence between fb/P 
from isolate 74 and fo/P from 75 and 76, all of which were STun- 
known48, is unlikely to have resulted from spontaneous mutation 
alone. The fo/P sequence encoded by isolate 74 was identical to 
that of isolates from four different serotypes, all of which were in 
the same clade. The presence of multiple STs in the same fo/P 
lineage is strongly suggestive of recombination. The fo/P sequence 
of ST989 of isolate 70 was identical to that of isolate 65. The high 
level of sequence divergence between folP from isolates belonging 
to ST989, STunknown22 and ST347 (7.2%, 5.6% and 7.3%, respect- 
ively) suggests that fo/P encoding DHPS with reduced affinity for 
sulfamethoxazole has spread via recombination. 

Our data show that, of the Malawian pneumococci that were 
trimethoprim resistant (MIC >32 mg/L), one-third possessed a 
previously unreported resistance genotype, a mutation at residue 
92 without Ile-100-Leu in DHFR, and suggest a potentially novel 
resistance mechanism. Previous studies report I100L as critical 
for the generation of trimethoprim resistance in pneumococci, 
while an additional mutation at residue 92 causes a further signifi- 
cant increase in MIC. 6 The substitution at residue 92 without the 
Ile-100-Leu mutation has been considered insufficient to confer 
trimethoprim resistance. 25 The majority of isolates (56%) pos- 
sessed Ile-100-Leu with an additional mutation at residue 92, 
and 10% contained Ile-100-Leu alone. Amongst those exhibiting 
the novel genotype there were no further amino acid substitutions 
in the trimethoprim-binding site. Amino acid substitutions outside 
the trimethoprim-binding site in DHFR were also found in the five 
trimethoprim-susceptible isolates, suggesting that they do not 
play a role in resistance. Trimethoprim-resistant pneumococci 
that do not possess the Ile-100-Leu substitution in DHFR have 
not been described previously. Further investigation is required to 
determine whether this novel pattern is unique to Malawi or has 
simply not yet been described in other geographical locations 
due to the small numbers of fo/A genes sequenced to date. 

Our data suggest that extensive co-trimoxazole use has biased 
the selection of pneumococci resistant to trimethoprim. The 



number of isolates with a mutation at position 92 but no 
Ile-100-Leu mutation and those encoding the Ile-100-Leu with 
the wild-type residue at position 92 decreased between 2004 
and 2008. In parallel there was a significant increase in those pos- 
sessing Ile-100-Leu with an additional mutation at residue 92. 
Over the same time period CPT increased substantially in the popu- 
lation. The additional mutation at position 92 has been shown to 
significantly reduce the affinity of DHFR for trimethoprim in relation 
to the Ile-100-Leu alone. 6 

Of the 20 STs represented more than once in the study popula- 
tion, only three shared an identical fo/A sequence and clustered to- 
gether in the phylogeny. ST did not predict how the remaining 17 
STs would cluster in the phylogeny, suggesting that recombination 
plays a major role in the spread of trimethoprim resistance. It is 
unclear why fo/A and fo/P are being disseminated by different 
mechanisms. 

There was an increase in the sequence divergence of fo/A 
amongst the isolates, matched by a decrease in dN/dS ratio from 
2004 to 2008. Neither change was significant, suggesting the se- 
quence divergence of fo/A from that of fo/A_R6 has remained con- 
stant. During the same time period the sequence divergence of fo/P 
from fo/P_R6 remained relatively stable, with minor drift. The years 
2002 and 2003 were excluded from the analysis because the 
sample contained fewer than three isolates from each of these 
years. The dN/dS ratio also remained stable from 2004 to 2008 
despite extensive co-trimoxazole use. Irrespective of isolation 
date, all resistant isolates displayed very high-level resistance 
to both trimethoprim and sulfamethoxazole. The sequence 
data therefore demonstrate that the evolution of both fo/A and 
fo/P is stable in the Malawian pneumococcal population. We pos- 
tulate that the strong selective pressure of extensive co- 
trimoxazole and sulfadoxine/pyrimethamine use over more 
than a decade means mutations additional to fo/A and fo/P will 
provide no additional fitness advantage, as the bacteria are 
now almost universally resistant. 

In conclusion, the employment of both an in vitro and an in silico 
analytical approach has enabled us to describe in detail the genetic 
basis of co-trimoxazole resistance in S. pneumoniae in Malawi. By 
creating robust phylogenies we have established that clonal dis- 
semination and recombination are contributing to the spread of 
co-trimoxazole resistance. Previously reported conventional 
testing for co-trimoxazole resistance in countries with high HIV 
seroprevalence has led to considerable debate as to why, with 
such high levels of pneumococcal resistance, clinical efficacy is 
maintained. 26,27 However, in the context of stabilizing evolution 
in silico, we suggest that the protective effects of co-trimoxazole 
in this vulnerable population may have changed since the relevant 
trials were conducted >10 years ago. 28,29 Our findings are likely to 
be generalizable to other countries where co-trimoxazole and 
sulfadoxine/pyrimethamine have been widely used. We therefore 
suggest that the efficacy of co-trimoxazole prophylaxis against 
pneumococcal infection in HIV/AIDS in sub-Saharan Africa 
should be re-evaluated. 
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